

Effects of Window 
Size and Shape on 
Accuracy of Subpixel 
Centroid Estimation 
of Target Images 


Sharon S. Welch 


(NASA-TP-3331) EFFECTS OF WINDOW 
SIZE AND SHAPE ON ACCURACY OF 
SUBPIXEL CENTROID ESTIMATION OF 
TARGET IMAGES (NASA) 37 p 


/Ay -3^ 
Sj'fi 


N9A-I5990 


Unc las 







NASA 

Technical 

Paper 

3331 


1993 


Effects of Window 
Size and Shape on 
Accuracy of Subpixel 
Centroid Estimation 
of Target Images 


Sharon S. Welch 
Langley Research Center 
Hampton, Virginia 


NASA 

National Aeronautics and 
Space Administration 

Office of Management 

Scientific and Technical 
Information Program 





Abstract 


A new algorithm, is presented for increasing the accuracy of subpixel 
centroid estimation of (nearly) point target images in cases where the 
signal-to-noise ratio is low and the signal amplitude and shape vary from 
frame to frame. In the algorithm , the centroid is calculated over a data 
window that is matched in width to the image distribution. Fourier 
analysis is used to explain the dependency of the centroid estimate 
on the size of the data window , and simulation and experimental 
results are presented which demonstrate the effects of window size for 
two different noise models. The effects of window shape have also 
been investigated for uniform and Gaussian- shaped windows. The 
new algorithm has been developed to improve the dynamic range of a 
close-range photogrammetric tracking system that provides feedback for 
control of a large gap magnetic suspension system (LGMSS). 


Introduction 

Centroid-estimation algorithms have long been 
used in digital imaging to locate target images to 
subpixel accuracies. Applications of centroid esti- 
mation include star tracking (ref. 1), point and edge 
detection for machine vision (ref. 2), close-range 
photogrammetry (ref. 3), and motion analysis. The 
effects of sampling and noise on the accuracy of the 
centroid estimate for point source images, images of 
extended sources, and edge detection have been ana- 
lyzed previously and documented by several authors 
(refs. 4 to 8). The systematic errors due to under- 
sampling that have been described for centroid esti- 
mation are common to all interpolation algorithms 
and have been analyzed from the point of view of 
performing image reconstruction (refs. 9 and 10). In 
these previous analyses, experimental approaches as 
well as analytical approaches based on Fourier tech- 
niques have been used to quantify the errors due to 
noise, quantization, and sample spacing. To date, 
the effect of window size on the accuracy of subpixel 
centroid estimation has been limited to a qualitative 
analysis derived from experiments that measured the 
error in centroid estimation as a function of different 
iV-point algorithms (ref. 4). 

In this paper, Fourier techniques are used to an- 
alyze the dependency of the systematic error on win- 
dow size. In addition, the effects of window size 
and shape on subpixel centroid-estimation accuracy 
in the presence of noise are studied. It is shown that 
there can be an advantage to using a shaped win- 
dow for centroid estimation of point target images for 
signals that vary in amplitude and width, provided 
the pixel-to-pixel noise is independent of signal am- 
plitude. A brief review of the effects of the optical 
point spread function (PSF), target size, and sam- 
ple spacing on systematic error is provided in order 


to compare and contrast these effects with those at- 
tributable to the data window. Quantization effects, 
however, are not addressed in this paper. 

In many applications involving centroid estima- 
tion, the signal shape and amplitude are either con- 
trollable or fixed. In these cases, the optimum sam- 
ple spacing and window size relative to the target 
image distribution are known a priori, and a correc- 
tion can be applied for systematic errors (ref. 7). In 
other applications w r here the noise is small or the im- 
age is averaged over several frames, a larger window^ 
(relative to the distribution) can be used in the cen- 
troid calculation. This eliminates systematic errors 
arising because of truncation of the signal. In the 
application discussed in this paper, centroid estima- 
tion is used to locate images of point targets along a 
linear charge-coupled device (CCD) detector, wrhere 
the signal-to-noise ratio is small and, because the tar- 
gets are moving, the (one-dimensional) images vary 
in size and amplitude. For these reasons, it is not 
possible to apply a correction for the systematic er- 
rors or to calculate the centroid using a large fixed 
data window. 

The application discussed herein is the optical 
measurement system (OMS) for the large gap mag- 
netic suspension system (LGMSS) (fig. 1). In the 
OMS, small infrared light-emitting diode (LED) tar- 
gets have been embedded in the top surface of a 
rigid cylindrically shaped element that contains a 
permanent magnet core. The element is magnetically 
levitated above a planar array of electromagnets. 
Sixteen linear CCD cameras arranged in pairs are 
located symmetrically about and above the levitated 
cylinder. A total of eight targets are located along 
the top surface of the cylinder (fig. 2), and the tar- 
gets are multiplexed in time for target identification. 
Triangulation techniques are used to determine the 


position and attitude of the levitated cylinder from 
the locations of the projected target images in the 
16 cameras. The position and attitude information 
is supplied to the electromagnet controller to stabi- 
lize levitation of the cylinder and to control motion 
in six degrees of freedom. 

The position and attitude of the cylinder are de- 
termined using weighted least squares. An estimate 
of the error in each computed centroid value is passed 
along with the centroid value and is used to estab- 
lish a weighting factor for the particular camera mea- 
surement. In order to achieve the required accuracy 
in the estimate of position and attitude, it is neces- 
sary to locate the centroids of at least 6 of the target 
images to Vis of a pixel in a minimum of 12 of the 
16 cameras. As the cylinder, and hence each tar- 
get, moves over the field of view, both the amplitude 
and the width of the target images vary (fig. 3). If 
the centroid location of a target image is determined 
with a fixed window size (which is the same for all 
16 cameras), and the window size is optimum for 
those image distributions falling in the midrange of 
possible values, then as the target images vary in am- 
plitude and width, the error in the centroid estimate 
grows for those images for which the distribution falls 
outside the midrange. Thus, with a fixed window 
size, the accuracy of the centroid estimate falls off be- 
cause of noise or systematic error, and the accuracy 
of the position and attitude estimate correspondingly 
decreases. 

In order to increase the dynamic range of the 
system, an algorithm has been developed to adjust 
the width of the centroid window as the light in- 
tensity distribution of the image varies. This algo- 
rithm provides the minimum error in centroid esti- 
mation over the maximum range of signal amplitude 
and shape. This paper analyzes the dependency of 
systematic and noise-induced errors on the size and 
shape of the data window. Experimental results are 
presented and compared with the results of numerical 
simulations. 

The following analysis is limited to one spatial di- 
mension. It is assumed that the PSF of the imaging 
optics can be approximated by a Gaussian function. 
Because a two-dimensional Gaussian function is sep- 
arable in x and y, the results of the one-dimensional 
analysis can readily be extended to two dimensions. 

Symbols and Abbreviations 

b half-width of Gaussian distribution 

used to simulate the optical point 
spread of the imaging optics 

comb(x/:r 5 ) sampling function 


/(*) 

light intensity distribution of one- 
dimensional target image 

m 

Fourier transform of f(x) 

m 

= mm 

f rw 

i 

■i 

* 

a; 

fe, 

II 

FswiO 

Fourier transform of fsw( x ) 

F SW r (0 

first derivative of with 

respect to £ 

fsw («) 

sampled and windowed version of 

m 

g{x- x 0 ) 

one-dimensional light intensity 
distribution of target located at x Q 

m) 

general function of £ 

N 

number of pixels corresponding to 
half the window width 

PSF(x) 

point spread function of imaging 
optics 

m 

Fourier transform of r(x) 

r(x) 

pixel response function 

no 

general function of £ 

w(0 

Fourier transform of w(x) 

w(x) 

window function 

X 

spatial variable 

X 

centroid of continuous one- 
dimensional light intensity 
distribution 

x c 

amount by which f(x) is shifted 
with respect to the sampling grid 

Xp 

pixel location of the peak signal 

x s 

sample spacing 

X SW 

centroid of fs\y{x) 

£c 

systematic error, x — xsw 

OVrH 

Abbreviations: 

variables of integration 

CCD 

charge-coupled device 

LED 

light-emitting diode 

rms 

root- mean-square 

SNR 

_ Peak signal 

Standard deviation of background signal 


A prime on a symbol denotes the first derivative. 
An asterisk used as an operational sign denotes con- 
volution of the quantities. 
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Effect of Window Size on Systematic 
Error 

The objective is to produce a minimum error 
estimate of the centroid of a one-dimensional light 
intensity distribution corresponding to a (nearly) 
point target image. The target image will have 
been sampled and only a portion of the sampled 
distribution (those points lying within the window) 
will have been used in the calculation. By definition, 
the centroid of the continuous one-dimensional light 
intensity distribution is given by 


set that is of infinite extent. In order to reproduce the 
real data set, which is of limited extent, the sampled 
function is multiplied by a window function u;(x) of 
finite width. If the sampled and windowed function 
is defined as fswi x )> then 

fsw = {f *r) (^-comb) w (4) 

where the asterisk denotes the convolution (pp. 284 
and 285 of ref. 11). Prom equations (3) and (4) it 
follows that the Fourier transform of fsw( x ) is 


_ _ I-OC x f( x ) dx 

X J^f(x)dx 


( 1 ) 


Fsw = F * comb * W 
= (F * W) * comb 


where x is the centroid, f(x) is the light intensity 
distribution of the image, and x is the position. This 
equation may be expressed in terms of F(£), the 
Fourier transform of /(x), as 


x = 


_—2niF(£) 


£=o 


( 2 ) 


where 


/ oo 

f(x) exp(— 27 rt£x) 

-OO 


dx 


( 3 ) 


The effect of using only a portion of the image 
distribution, that is, of placing a window about the 
distribution, is to spread the energy in the spectrum 
of F(£) into higher frequencies. This effect is shown 
in figure 4. When the image distribution is both 
windowed and sampled, as it is in practice, the width 
of the window relative to the width of f(x) and the 
sample spacing is important. 

In a CCD array, an image is sampled at discrete 
intervals with an array of detectors or pixels. Each 
pixel has a response that is spatially distributed over 
the width of the pixel. If f(x) is the function that 
represents the continuous light intensity distribution 
of the image, then sampling with a CCD array can 
be modelled mathematically as a convolution of f(x) 
with a function r(x), which represents the spatial 
response of each pixel, followed by a multiplication 
of the resulting product with the sampling function 
comb(x/x s ), where 


— comb(x/x 5 ) 
x s 


oo 

^ 6(x — nx s ) 

n— —oo 


and x s is the sample spacing of the pixels. This con- 
volution and multiplication results in a sampled data 


where 


and 

Thus, 


comb = 


E « U 


n 

Xs 


F = FR 


F sw(0 = E F nw ( f ~ ~r) 

v X S / 


( 5 ) 


where _ 

Frw = F*W (6) 

Equation (5) shows that the spectrum of the 
sampled and windowed function is the sum of the 
individual spectra Fp^y repeated at intervals of 
l/x s . Figure 5 shows a representative sampled and 
windowed function fsw( x ) and the corresponding 
Fourier transform Fsw/(£)* As for the continuous 
distribution case, placing a window about the image 
distribution can broaden Fji\y. When the image is 
both windowed and sampled this spread can result 
in energy spillover from neighboring transforms. If 
the window width is too small relative to the width 
of the signal and the sample spacing, then system- 
atic error is introduced into the centroid estimate as 
energy from higher orders spills into Fr\y( 0). 

Before the specific form of the error introduced 
by the data window is discussed, the effects of the 
optical PSF, pixel response function, and target size 
are briefly reviewed in order to distinguish between 
those errors that are due to undersampling and those 
errors that are due to the finite extent of the data. 

A Review of Effects of Optical PSF, Pixel 

Response Function, and Target Size 

As has been shown previously (refs. 4, 7, 9, 
and 10), the form of the subpixel error in centroid 
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estimation and image reconstruction due to under- 
sampling is a sinusoidal function of the position of 
the “true” centroid or image location with respect to 
the sample grid (fig. 6). The magnitude of this error 
is a function of the widths of the optical PSF, pixel 
response function, and target. 

In a digital imaging system, the optical PSF 
describes the degree to which the target image is 
blurred because of the limits of diffraction in the 
imaging optics. Provided the imaging optics are shift 
invariant, then the light intensity distribution of the 
image f(x — x Q ) is the convolution of g(x — x a ) with 
the PSF of the imaging optics, where g(x — x 0 ) is 
the light intensity distribution of the true target and 
x 0 is the location of the centroid of the distribution. 
(Sec pp. 335 and 336 of ref. 11.) 

Because the optics tend to blur, or smear out 
the image of the target, the PSF has an inverse 
effect on the Fourier transform of the image. That 
is, the Fourier transform of the blurred image is 
narrower than the Fourier transform of the true 
image of the target. For this reason, as shown in 
references 4 and 6 to 8, the systematic error in 
centroid estimation, or image reconstruction, due to 
undersampling decreases as the width of the PSF of 
the imaging system increases. 

Similarly, the pixel response function r(x) spreads 
the sampled image distribution. This spreading is 
shown in equation (4), where r(x) is convolved with 
f(x ), or equivalently in equation (6), where F(£) 
is multiplied by 7?(£). For the linear CCD array 
detector that was used to generate the experimental 


results discussed in this paper, the width of the pixel 
response function is approximately equal to Q.8x s . 

The size of the target also affects the amplitude 
of the undersampling error. Typically, the larger the 
target relative to the sample spacing, the smaller 
the systematic error. However, the amplitude of the 
undersampling error does not decrease monotonically 
with increasing target width. Rather, for a fixed pixel 
response function, if the amplitude of the error is 
plotted as a function of target size, then the resulting 
graph shows this amplitude to be modulated period- 
ically with target size (refs. 5 and 8). This period- 
icity can be understood by looking at the graphs of 
figure 7, wherein simulated functions f(x) and F(£) 
are plotted for two different target widths, x s and 4x s 
(uniform amplitude). The PSF (fig. 7 (a)) is assumed 
to be of the form 

PSF(x)=exp -*(^) (7) 

As is evident in the graphs of figures 7(b) to 7(d), if 
the target is uniform in intensity and has sharp dis- 
continuities at the edges, then the Fourier transform 
of the true target image looks like a sine function 
(fig. 7(d)). As the width of the target changes for a 
fixed PSF and pixel response function, the side lobes 
move relative to zero frequency. For some target sizes 
and PSF’s, F(£) and F'(£) are zero; for others, F'(£) 
is between zero and a peak of the side lobes. The 
result of this is that the amplitude of the sinusoidal 
error as a function of target size is modulated in a 
periodic fashion that is determined by the ratio of 
the width of the target to the sample spacing. 


r 


t 


Centroid-Estimation Error for Sampled and Windowed Function 

From equation (2), the estimate of the centroid of the sampled and windowed function f$w( x ) ls 


X SW = 


4(0) 

—2niFgw{fy 


( 8 ) 


With the systematic error e c defined as the difference between x and xsw ■ then from equations (2) and (8) 


i l>(0) EswM 

e c -x xjw 2 ^[ F(0) f sw(0 )_ 


( 9 ) 
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Substituting equation (5) for yields 


e c 


1 

2ni 


F'(Z) 

m 


oo 

(d/di) E F RW ((, - (n/*s)) 

n=— oo 

E f rw (£ - ( n / x s)) 

n=—oo 


e=o 


( 10 ) 


Using the approach introduced in reference 7, and assuming that the distribution of the target image is an 
even function about a point that is shifted an amount x c with respect to the sampling grid, results in 


/ 0*0 = fe{x-x c ) 


and 

F(£) = exp(-2iriXcO F e{(>) 

Carrying this analysis further, the effect of the window on the Fourier transform F(£) is, as shown by 
equation (6), a convolution of the Fourier transform of the window function VF(£) with F(£), or 

Frw =F*W= |exp( — 2irix c fi ) F e J *W = J cxp(-2irix c r,)F c (r,)W (£ - r,) dr] (11) 

To see the effect of the window on the error e c , we return now to equation (10) and solve for F^jy(£) using 
equation (11) and the following relationship: 

/ oo roc 

T( V )H(Z - r,) dr, = / T(i- r,)H(r,) dr, 

-OO J —oo 

This yields 

/ oo _ 

exp [-2? rix c (£ - rj)] F e (£ - r,)W(r,) dr, 

-OO 


and 


d roc _ 

w(0 = -77 / ex P [-27 rix c {£, - r,)} F e (£ - r,)W{r,) dr, 

J— OO 


= —2nix { 


T 


/ exp [—2rrix c (£ - 77)] F e (£ - r,)W(r,) dr, 
J —00 

/ oo _ 

exp[-27rix c (£ - 7?)]Fg(^ - 77)1^(77) dr, 

-OO 


Plugging these values for F/^(£) and F^y(^) into equation (10), setting £ — 0, and cancelling like terms in 
the numerator and denominator leaves 


E (f^ 0 00 exp{-2nix c [-r,-(n/xs)]}F/,[-r}-{n/x s )]W(r])dr,^ 
1 n=— oo ' 

2m E (/^ 0 exp{- 27 rix c [- 7 ? - (n/x s )}} F e [-r, - (n/x a )]W(r,) dqj 

71= — OO ' 


( 12 ) 


If we assume that the window is an odd number of pixels, then VF(£) is an even function. Furthermore, if 
one makes the assumption that the pixel response function is symmetric about the center of the pixel, then 
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R{€) is also an even function. Rearranging terms in the numerator and denominator of equation (12) and 
realizing that F f e is an odd function yields 


(- 


€ r = 


27 r 


/-°°oc s M 27r x c ri)Fe(r))W e (ri) dr] + 2 E {/^oSin(27rx c M)F'(/i)W e [(n/x s ) - n] d/x} 
I^oo cos(2nx c Tj)F e (r))W c (i]) dr) + 2 E {/^o cos(27rx c /x)F e (^)W e [(n/x s ) - /x] d/x) 


(13) 


where 


and 


/oo _ 

cos(27rj; c /i)Fg(//)W e (/i) dfi = 0 

-OC 

/ oo _ 

sin(27 Tx c ^)F e (fi)W e (fi) d^i — 0 

-oo 

because the integral of an odd function times an even function evaluated over even limits is zero. The first 
integral in the denominator of equation (13) is generally much larger than the second integral, and therefore 


/. 


2tt 


J™oc sin(27rx c 77)F e '(^)W e (7y) dr] + 2 £ |/^ 0 sin(27rx c /x)F'(/x)W r e [(n/x s ) - /x] d/x) 


/f^ 0 cos(27rx c j?)F e (7?)W e (r7) dr? 
When the window is large, W e (£) — ► <5(£) and equation (14) reduces to 


v 


(14) 




1 


oo _ 

E sin(27 Tx c njx s )F' e {n/x s ) 

71=1 


W) 


Because of the frequency cutoff of the optics, the n — 1 term dominates. Thus, the form of the systematic 
error in this case is sinusoidal with an amplitude proportional to F f e {\jx s ). 

However, when the width of the window is approximately the same as the width of the image distribution, 
the convolution integrals of equation (14) spread each of the individual Frw spectra and terms higher than 
the first become significant. As a result, the form of the subpixel error changes from a pure sinusoid to one 
period of a sawtooth plot (fig. 8). The sawtooth form corresponds to the systematic error that arises because of 
truncation of the signal. The amplitude of this error is now a function of the width of the window in addition 
to the widths of the pixel response function, the PSF, and the target. 

A simulation was constructed to study the relationship between the shape and width of the window and 
the systematic error in centroid estimation. In all the simulations discussed herein, the target was assumed to 
be a point source, the PSF was assumed to be that given in equation (7), and the pixel response function was 
assumed to be uniform and of width x s . Two different-shaped windows were simulated: (1) a uniform window, 
defined by 

{ > f 1 (x v — Nx s < x < x v + Nx s for TV = 1,2,...) 

= { 0 (Otherwise) 

and (2) a Gaussian-shaped window defined by 


( \ _ / cxp[-7r(x - Xp)/(2Nx s ) 2 } (x p — 2Nx s < x < x p + 2Nx s for N = 1, 2, . . .) 

W{X) ~ 1 0 (Otherwise) 

where x p is the pixel location of the peak signal. The functions w(x) and corresponding Fourier transform 
W(£) for the uniform and Gaussian-shaped windows are shown for comparison purposes in figure 9. Note that 
the Fourier transforms for both windows have a zero at approximately the same location. The width of the 
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Gaussian-shaped window has been chosen to provide a similar shape in frequency space to that of the uniform 
window. In addition, the error due to undersampling has been made negligibly small. Provided b > 1.25x 5) 
which was true for the cases th^,t were simulated, the maximum error due to undersampling is less than 0.002 
times the pixel spacing. 

In figure 10, the simulation results show the predicted peak amplitude of the subpixel systematic error as 
a function of the ratio N/b for the two different- shaped windows. These results show that using a Gaussian- 
shaped window reduces the systematic error. The systematic error due to windowing goes to zero with the 
Gaussian-shaped window for values of N/b greater than approximately 1.0 and with the uniform window for 
values of N/b greater than approximately 2.5. 


Effects of Window Size and Shape on 
Noise-Induced Error 

In order to understand the effects of window size 
and shape on centroid-estimation error in the pres- 
ence of noise, a simulation was constructed. Two 
different noise models were simulated. In one, the 
noise was modelled as a mean zero Gaussian ran- 
dom variable, independent of signal amplitude and 
independent for each pixel along the array. In the 
other, the noise was modelled as being proportional 
to the square root of the signal amplitude, where the 
proportionality factor was modelled as a mean zero 
Gaussian random variable. The reasons for selecting 
these two noise models are the following. Often in 
imaging applications the target signal is small rela- 
tive to the background light, and the noise is con- 
sidered to be background limited photon shot noise 
(BLIP). This case is simulated by the noise model 
in which the noise is independent of the signal and 
independent pixel to pixel. The other case closely 
approximates the characteristics of the sensor noise 
in the OMS. 

Equation (7) was used to mathematically model 
the PSF of the optics. The source was modelled as a 
delta function with a peak signal of 1.0. A constant 
background signal was added to the target signal 
prior to addition of the noise. When the centroid was 
estimated, the values at the endpoints of the window 
were averaged, and this average value was subtracted 
from each signal value. This process simulated the 
threshold technique that is used on the OMS. 

Simulations were run to determine the subpixel 
bias error and the standard deviation of the signal 
at any single subpixel position. To calculate the bias 
error, an average centroid estimate was calculated at 
each subpixel position and subtracted from the true 
centroid location. The average centroid estimate was 
taken to be the mean of 25 samples. The standard 
deviation of the centroid estimate at each subpixel 
location was also calculated. 


The results of this simulation analysis are shown 
in figures 11 to 13. In figures 11(a) to 11(d), the av- 
erage error over a single pixel is plotted as a function 
of x c for representative samples of each combination 
of noise model and window shape. In figures 12(a) 
and 12(b), the rms value of the bias error over a sin- 
gle pixel is plotted as a function of the ratio N/b for 
the uniform window and the Gaussian-shaped win- 
dow with the pixel-to-pixel noise independent. The 
rms error over the range —0.5 < x c < 0.5 is defined 
as 


rms = 



P 

^"^(True centroid — Predicted centroid) 2 

71=1 


where P is the total number of subpixel positions. 
In figures 12(c) and 12(d), similar plots are shown 
for the noise proportional to the square root of the 
signal amplitude. The different points correspond to 
different values of SNR in figures 12(a) and 12(b), 
where the SNR is defined as 

Peak signal 

Standard deviation of background signal 

and to different standard deviations of the propor- 
tionality factor in figures 12(c) and 12(d). Fig- 
ure 13(a) presents the standard deviations of the 
centroid estimates at single subpixel locations as a 
function of the SNR for the uniform and Gaussian- 
shaped windows and noise independent of signal am- 
plitude. In figure 13(b), similar graphs are shown for 
the uniform and Gaussian-shaped windows and noise 
proportional to the square root of signal amplitude. 

Four things become apparent from the simula- 
tion results, the first of which has been known for 
some time: (1) in the presence of noise, the standard 
deviation of the subpixel centroid estimate increases 
with increasing window size, no matter what the win- 
dow shape or noise model; (2) the optimum window 
size, that is, the window size that minimizes both the 
bias error and the standard deviation of the centroid 
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estimate, is a function of the SNR and the noise 
model; (3) the optimum window shape is a function 
of the noise model; and (4) the subpixel bias error 
tends toward the sawtooth shape associated with the 
systematic truncation error when the noise is pro- 
portional to the signal. For both noise models, the 
rms value of the bias error increases at a slower rate 
for the Gaussian-shaped window than for the uni- 
form window. However, the magnitude of this error 
is greater for the Gaussian-shaped window than for 
the uniform window when, the noise is proportional 
to the square root of the signal amplitude; the error 
magnitude for the Gaussian-shaped window is less 
than that for the uniform window when the pixel-to~ 
pixel noise is independent and mean zero Gaussian. 
It follows from this that in order to determine an 
optimum window size (and shape) for point target 
tracking, the SNR and noise process must be known. 

Experimental Results 

Experiments were run to verify the results pre- 
dicted by simulation. A single camera was mounted 
directly over an LED target. The target was moved 
over ±0.05 in. in the object plane in 0.0025-in. steps 
using a computer-controlled linear stage and imaged 
with the camera (fig. 14). The displacement of the 
stage was measured using a laser interferometer to 
approximately 0.0001 in. Since the magnification fac- 
tor for the camera was approximately 0.1, it was pos- 
sible to resolve the displacement of the target in the 
image plane to better than V 20 of a pixel. The cen- 
troid of the projected target image was calculated 
as the average of the centroid estimates obtained for 
100 target images acquired at each location of the 
stage. A straight line was fit through the average cen- 
troid estimates as a function of the measured stage 
position. The bias error was then determined to be 
the difference between the straight-line fit and the 
average centroid estimate at each location. 

Figures 15(a) and 15(b) show the bias errors for 
a uniform window and a Gaussian-shaped window. 
In figure 15(a) (for the uniform window), the larger 
amplitude error corresponds to N = 4 and the 
smaller amplitude error corresponds to N — 6. In 
figure 15(b) (for the Gaussian-shaped window), the 
larger amplitude error corresponds to N = 5 and 
the smaller amplitude error corresponds to N = 15. 
The image distribution for which the centroid was 
estimated is that shown in figure 16 and was the same 
for all N. Assuming a point target, the half-width 
of the distribution shown in figure 16 corresponds to 
b & 6.0. The jitter that is apparent in the data is 
the result of the noise model and the thresholding 
technique discussed previously. 


In order to determine the relationship between 
the signal SNR, window size and shape, and accu- 
racy, the experiment described above was repeated 
for different sensor integration times. Figure 17 
presents the rms values of the bias error (over 9 cy- 
cles) for different SNR’s and window widths. Fig- 
ure 17(a) shows the results obtained with a uniform 
window, and figure 17(b) shows the results obtained 
with a Gaussian-shaped window. 

In addition to the bias error, the standard de- 
viation of the centroid estimate at a single position 
of the stage was measured. These results are shown 
in figure 18. As expected, the standard deviation of 
the centroid estimate increases with increasing win- 
dow size. The increase in the standard deviation 
again at small window sizes is believed to result from 
thresholding. As the window size decreases to the 
point where the signal is truncated, those pixels at ei- 
ther edge of the window have greater noise associated 
with them. When the threshold is set as the average 
of the signals from these two pixels and subtracted 
from the signal at each pixel within the window, the 
random error in the estimate of the centroid is in- 
creased. Therefore, the variation from sample to 
sample increases. 

From the above results, the uniform window 
shape is optimum for centroid estimation in the op- 
tical tracking system. Furthermore, N/b « 1.5 is the 
optimum ratio of window width to width of the target 
image distribution for SNR’s from about 50 to 100. 
This ratio of window size to distribution width pro- 
vides the minimum total error, where the total error 
is the sum of the rms value of the bias error and the 
standard deviation of the centroid estimate. 

Centroid-Estimation Algorithm 

A flowchart of the centroid-estimation algorithm 
is shown in figure 19. The subpixel centroid estimate 
is calculated in the following way. During tracking, a 
search is made for the peak pixel location of each tar- 
get prior to computation of the centroid. The pixel 
location of peak light intensity is stored and used 
to set up a window of 81 pixels centered about the 
peak. Each pixel has a slightly different response, 
that is, a slightly different gain and zero point. Prior 
to computation of the centroid, the values of light 
intensity are corrected for pixel nonuniform respon- 
sivity and the background light level is subtracted 
for pixels falling within the window. A second coarse 
search for the peak intensity is then performed over 
the corrected light intensity values. 

Before the centroid is estimated, the width of the 
distribution is determined by computing the second 
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moment of the target image distribution. The second 
moment is calculated as 



where x is the centroid of the distribution calculated 
over ±15 pixels centered about the peak, I in) is the 
digital value of the light intensity for pixel n, and x p 
is the pixel location of the peak signal. 

The value of N used in the centroid calculation 
is equal to the nearest integer value of the second 
moment a plus one pixel. The subpixel centroid 
estimate is calculated as 

Xp+N 

E nI ( n ) 

n~Xp — N 
Xsw = x p +N 

E I in) 

n—Xp—N 


Figure 20 shows the subpixel bias error from the use 
of the centroid-estimation algorithm that adjusts the 
window size; the window shape is uniform. Figure 21 
shows the standard deviations of the centroid esti- 
mate for a uniform window of variable width and 
a uniform window with a fixed width of N = 14. 
The value of N = 14 was chosen as the size required 
to minimize the systematic error for the maximum 
signal obtained during tracking with the OMS. The 
maximum signal has an SNR of approximately 200. 
Adjusting the window size to match the width of the 
distribution minimizes the systematic error as well as 
the standard deviation of the centroid estimate over 
a wider range of SNR's than does using a window of 
fixed size. 


Conclusions 

The effects of window size and shape on the ac- 
curacy of subpixel centroid estimation have been 
presented. Two different noise models and win- 
dow shapes have been studied. The shapes include 
a uniform window and a Gaussian- shaped window. 
The noise models studied include random mean zero 
Gaussian, independent pixel to pixel and indepen- 
dent of signal amplitude, and noise proportional to 
the square root of the signal amplitude. 


Fourier analysis has been used to determine the 
form and magnitude of the systematic error due 
to windowing. It has been shown that the form 
of the subpixel error due to windowing resembles 
one period of a sawtooth compared with the pure 
sinusoidal waveform that is obtained when the signal 
is undersampled. Furthermore, it has been shown 
that the magnitude of this systematic error is smaller 
for a Gaussian-shaped window^ than for a uniform 
window. 

Simulations have been run to explore the sensitiv- 
ity of the centroid-estimation algorithm to window' 
size and shape in the presence of noise and to deter- 
mine the window shape and size that minimize the 
error. Experiments have been conducted to verify 
the behavior predicted by simulation. The results of 
the simulations and experiments revealed the follow- 
ing: (1) for noise that is proportional to the square 
root of the signal amplitude, the optimum wundow 
shape is uniform, and for noise that is independent 
of signal amplitude, with a low ratio of peak signal to 
standard deviation of background signal (SNR), the 
optimum window shape is Gaussian; and (2) match- 
ing the size of the window' to the width of the target 
image distribution improves the accuracy of the cen- 
troid estimate for both noise models and both win- 
dow shapes. The optimum ratio of window^ width 
to width of the target distribution is approximately 
1.5 for the uniform window' and noise proportional to 
signal amplitude. 

The results of the analysis have been used to de- 
velop a new centroid-estimation algorithm that in- 
creases the accuracy of subpixel centroid estimation 
of (nearly) point target images w'hen the noise is pro- 
portional to signal amplitude and the signal ampli- 
tude and shape vary from frame to frame. In the al- 
gorithm, the width of the data window' is matched to 
the estimated width of the image distribution. Cal- 
culating the centroid over a window that is matched 
in size to the width of the distribution yields a sub- 
pixel centroid estimate with smaller total error over 
a wider range of SNR’s than that from a calculation 
with a window of fixed size. This improvement in 
centroid estimation has been developed for a point 
target tracking system in order to increase the dy- 
namic range of the system. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
June 22, 1993 
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Figure 1. Optical measurement system (OMS) of large gap magnetic suspension system (LGMSS). 




Rechargeable batteries 

Figure 2. Levitated cylinder with locations of eight LED targets. 
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Digital signal amplitude 


i 


300 



Pixels along array 


Figure 3. Digital signals (expanded) corresponding to light intensity distributions for images of targets 1, 3, 
5 ? and 7 as viewed from sensor 7 in OMS. All images have been referenced to same peak pixel location to 
better illustrate range of sizes and shapes of sampled image distributions that can be generated by different 
rgete 


Spatial frequency 

(a) Continuous distribution f{x). 


x 



(b) Fourier transform of /(#) multiplied by Fourier transform of window function w(x), uniform window placed 
symmetrically about f(x). 


Figure 4. Effect of windowing. 



Spatial frequency 

(b) Fourier transform of signal. 

Figure 5. Representative sampled and windowed signal corresponding to point target image and Fourier 
transform of this signal. 





Figure 6. Error in centroid estimate as function of x c generated from simulated data. Target was assumed to 
be point source, PSF was assumed to be Gaussian with widths ranging from 1.0 to 1.5 times sample spacing 
x s . and pixel response function was assumed to be uniform and equal in width to sample spacing. 






Systematic error, pixels 





Amplitude Amplitude 




(b) W(Z). 

Figure 9. Representative functions w(x) and VF(£) for Gaussian-shaped window and uniform window. N = 5 






(a) Gaussian mean zero noise and uniform window 
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(c) Noise proportional to square root of signal 


amplitude and uniform window. 
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(a) Gaussian mean zero noise. 
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(b) Noise proportional to square root of signal amplitude. 
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Figure 13. Predicted standard deviations of centroid estimates at any single subpixel position. 
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Figure 14. 
stages. 


Experimental setup showing camera mounted overhead of diode and computer-controlled linear 
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(a) Uniform window. 



(b) Gaussian-shaped window. 


Figure 17. Experimental rms of bias error as function of window size. Note the difference in vertical scale 
between the two plots. 




















Bias 





o Fixed width 
° Variable width 






REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503 


1. AGENCY USE ONLYfleave blank) 2. REPORT DATE 

October 1993 


4. TITLE AND SUBTITLE 

Effects of Window Size and Shape on Accuracy of Subpixel 
Centroid Estimation of Target Images 


REPORT TYPE AND DATES COVERED 

Technical Paper 


5. FUNDING NUMBERS 

WU 590-14-11-02 


6. AUTHOR(S) 

Sharon S. Welch 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Langley Research Center 
Hampton, VA 23681-0001 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 

L-17113 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA TP-3331 



12a. DISTRIBUTION/AVAILABILITY STATEMENT 


12b. DISTRIBUTION CODE 


Unclassified - Unlimited 


Subject Category 35 


13. ABSTRACT ( Maximum 200 words) 

A new algorithm is presented for increasing the accuracy of subpixel centroid estimation of (nearly) point 
target images in cases where the signal-to-noise ratio is low and the signal amplitude and shape vary from 
frame to frame. In the algorithm, the centroid is calculated over a data window that is matched in width 
to the image distribution. Fourier analysis is used to explain the dependency of the centroid estimate on the 
size of the data window, and simulation and experimental results are presented which demonstrate the effects 
of window size for two different noise models. The effects of window shape have also been investigated for 
uniform and Gaussian-shaped windows. The new algorithm has been developed to improve the dynamic range 
of a close-range photogram metric tracking system that provides feedback for control of a large gap magnetic 
suspension system (LGMSS). 


14. SUBJECT TERMS 

Centroid estimation; Target tracking; Charge-coupled device (CCD) sensors 


15. NUMBER OF PAGES 

36 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 




18. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 20. LIMITATION 
OF THIS PAGE OF ABSTRACT OF ABSTRACT 

Unclassified Unclassified 


Standard Form 298(Rev. 2 

Prescribed by ANSI Std Z39-18 
298-102 























